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While different ways to realize spin-orbit coupling in Bose-Einstein condensates exist, not all are 
currently experimentally implementable. Here we present a detailed study of gap solitons in a Bose- 
Einstein condensate with experimentally realizable spin-orbit coupling and discuss two cases relating 
to a spin-dependent parity symmetry. In the parity symmetric case, two families of fundamental 
gap solitons in second linear energy gap are demonstrated with opposite sign of the parity, with 
one family having single humped densities and the other double humped ones. In the case of 
broken parity symmetry, the fundamental modes manifest spin-polarization. Both families possess 
an opposite sign of the spin-polarization. 
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I. INTRODUCTION 

The interplay between spatial periodicity and non¬ 
linear dynamical evolution harbours a large number of 
interesting physical many body effects. A particularly 
outstanding one is the existence of unique soliton solu¬ 
tions that reside in the linear energy gaps between the 
Bloch bands [1, 2]. These so-called gap solitons are char¬ 
acterized by several remarkable features: independent of 
the repulsive or attractive nature of the nonlinearity, gap- 
solitons are always bright solitons [3] and even though 
they are spatial localized, they are intrinsically related 
to the extended Bloch waves [4-7]. The general prop¬ 
erties of gap solitons have been intensively studied in 
many different nonlinear periodic systems, such as opti¬ 
cal waveguide arrays and photonic lattices [1], semicon¬ 
ductor microcavities [8, 9] and also atomic Bose-Einstein 
condensates (BECs) in optical lattices [10]. The latter 
systems allows very good control over all available phys¬ 
ical parameters and matter-wave gap solitons [11] have 
been observed experimentally for the first time in 2004 
[ 12 ]. 

Another recent area of interest related to atomic BECs 
is the behaviour of the atoms in the presence of spin-orbit 
coupling [13-17]. Eor neutral atoms this corresponds to a 
coupling between a pseudo-spin and the motional degrees 
of freedom and this interplay has in recent years shown 
to lead to intriguing new phases [18-24] and to modified 
excitations and dynamics [22, 25-29] (see [30-34] for re¬ 
views). Similarly, solitonic solutions in the presence of 
spin-orbit coupling are predicted to show many interest¬ 
ing properties [35-43]. However, as one of the prominent 
features of spin-orbit-coupled system is the presence of a 
broken Galilean invariance [44] , the conventional method 
of finding moving solitons using stationary ones is not ap¬ 
plicable. It is therefore a nontrivial task to find moving 
soliton solutions [38, 43]. 

The combination of spin-orbit coupling and optical lat¬ 
tices in BEC experiments is particularly intriguing. Spin- 
orbit coupling dramatically changes the Bloch spectrum 
and, under suitable parameters, the lowest Bloch bands 


can become flat [45, 46]. The collective excitation pos¬ 
sesses novel roton-like structures [47]. These new spec¬ 
tra are known to modify the standard Mott-superfluid 
transition [48-50] and allow for the existence of many 
exotic magnetic states [51, 52]. Eurthermore, nontrivial 
dynamics, such as unconventional Bloch oscillation [53] 
and anomalous Anderson localization [54-56], are pre¬ 
dicted to appear. 

Recently, Kartashov et al. studied the properties of 
gap solitons in a spin-orbit-coupled BEC in one and two 
dimensional spin-dependent optical lattices and demon¬ 
strated that gap solitons can be categorized according 
to their spin-dependent parity, time-reversal and trans¬ 
lational symmetries [57, 58] . These studies provide a use¬ 
ful guideline for finding gap soliton solutions. However, 
the discussed systems have not yet been experimentally 
realized. 

So far, only two experimental approaches to prepare 
spin-orbit-coupled BECs in lattice potentials exist. The 
first is a combination of Raman lasers and a radio¬ 
frequency magnetic field [59], where both the Raman 
beams and the radio-frequency field couple two ener¬ 
getic ground states of the atoms. The momentum im¬ 
parted to the atoms from the Raman beams generates 
artihcial spin-orbit coupling, and meanwhile dresses the 
radio-frequency coupling to create spin-dependent lattice 
potentials. Eor this, however, the frequency difference be¬ 
tween the Raman beams has to be exactly the same as 
the frequency of the radio-frequency field, which is a sig¬ 
nificant limitation for the parameters space that can be 
explored. The second approach is to adiabatically load a 
Raman coupled BEC into a conventional optical lattice 
that is formed by two counter-propagating laser beams 
[60]. A frequency difference between the lattice beams 
creates a moving lattice and by tuning the frequency dif¬ 
ference, the BEC can be loaded into any Bloch band in 
the moving frame. 

Motivated by the experimental developments reported 
in Ref. [60] and theoretical studies [57, 58], we present 
in this work a systematical study of the properties of 
gap solitons in an experimentally realisable spin-orbit- 
coupled BEC in a one dimensional optical lattice. The 
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properties of gap solitons existing in the first linear gap 
for the repulsive interactions were first studies in Ref. [57] 
and here we mostly focus on the two families of funda¬ 
mental modes in the second gap, which have a number 
of interesting properties. If the detuning between the 
Raman beams and energy states of the atoms is zero, 
the fundamental modes have well-defined parity symme¬ 
try and different families of the fundamental modes take 
the opposite sign of the parity. Furthermore, the density 
distributions of the second family has a double humped 
structure, which disappears in the presence of small Ra¬ 
man coupling, when the density distribution of the sec¬ 
ond family becomes similar to that of the first family. 
For non-zero detuning, the Hamiltonian does not have 
the parity symmetry and we find that the spin composi¬ 
tions of different fundamental mode families differ signif¬ 
icantly. This provides a possible experimental signature 
to distinguish them. 

This manuscript is organized as follows. In Sec. II we 
introduce the description of our system in terms of di¬ 
mensionless Gross-Pitaevskii (GP) equations. After that, 
in Sec. Ill, we present the linear Bloch spectrum and 
Bloch waves, from which the position of the linear en¬ 
ergy gaps are identified. The parity symmetry of the 
Bloch waves is illustrated. In Sec. IV we demonstrate 
in detail the properties of gap solitons considering the 
case of zero detuning. The case of non-zero detuning is 
analyzed in Sec. V and we conclude in Sec. VI. 


II. MODEL 


Our study is fully based on the very recent experimen¬ 
tal realization of a spin-orbit-coupled BEG by Hamner et 
al. [60]. To create spin-orbit coupling, they select two 
energy levels of ®^Rb atoms, which act the pseudo-spin- 
1/2, and Raman couple them with a coupling strength 
H, called the Rabi frequency, and with a detuning given 
by 5. This Raman coupling creates the spin-orbit cou¬ 
pling through the momentum transfer between the Ra¬ 
man beams and the atoms and the Hamiltonian that re¬ 
sults from this coupling is given by 

rr PI , I I ^ 

tisoc = -1- ^PzCfz + TTCTz -I- —Ox- (1) 

2 m 2 2 

Here the {(Ji} are the Pauli matrices, m is the atom mass, 
and Pz is the momentum of the atoms along the Ra¬ 
man beams, which counter propagate along the 2 direc¬ 
tion. Spin-orbit coupling happens in the same direction 
and the coefficient 7 = is proportional to the 

wavenumber of the Raman beams /cRam- 

The spin-orbit-coupled BEG is then loaded into a one 
dimensional optical lattice by adiabatically ramping up 
the lattice beams, which are optically aligned with the 
Raman beams. The mean-field description of the spin- 
orbit-coupled lattice BEG is then given by the coupled 


GP equations, 

— [T^soc “f V sin (A^lat-^)] d/ -t- 

( 2 ) 


where the nonlinear terms H^on originate from the atomic 
interactions. They are explicitly given by 

an [TilVoiz 1 «' 2 |^ 0 

0 ai 2 l«'ilVa 22 l«' 2 |^ 

(3) 

and the spinor wave-function is dl = (vki,'I' 2 )^. Here, 
/ciat is the wavenumber of lattice beams, v is lattice depth 
which can be tuned by changing the intensity of the lat¬ 
tice beams, and ujx and ojy are the harmonic trap fre¬ 
quency in the transversal directions. We study an elon¬ 
gated BEG, so that the weak harmonic trap along the z 
direction can be neglected, and we assume that the dy¬ 
namics in the transverse directions is completely frozen 
out due to large transversal trapping frequencies. The 
®^Rb atoms are characterized by repulsive interactions 
and the relevant s-wave scattering lengths are quanti¬ 
fied by 011,022 and 012 . In the experiment by Hamner 
et al. [60] the difference between the s-wave scattering 
lengths is very small and in the following we will assume 
Oil = 022 = 0-12 = a for simplicity. 

To find gap solitons numerically, we start from the di¬ 
mensionless GP equations 




I 52 .5 <5 n ' 


$ 


- ^ cos( 2 z)$ -h (l$il^ -k 1 ^' 2 |^)^', 


(4) 


where the units of energy, time and length are 
2 Eiat = Ii/ 2 Eiat and I/fciat respectively. 

The dimensionless wave-function is given by $ = 

T so that the number of atoms N turns 

out as 


N = Noydz(l$il" + l$ 2 |"), (5) 

with 

Nq — Plat / (^\/ OJx^y ^latC^) ■ (6) 

For the detailed numerical calculations presented below 
we use the experimental parameters given in [60], fciat = 
27r/V2A lat with Alat — 154071771 and AlRam — 27r/'\/2ARam 
with ARam = 784nTO, which leads to a dimensionless spin- 
orbit coupling parameter 7 = fcRam/Mat = 1.96. Through 
out this work we set v = 2 corresponding to 4Eiat. The 
detuning S and the Rabi frequency H are free parameters 
and in the experiment they can be tuned by the chang¬ 
ing of the frequency difference and intensity of Raman 
beams, respectively. 

Gap solitons are stationary solutions of Eq. (4) and 
have the form ^{z, t) = <j){z) exp{—ifj,t), with p being the 
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FIG. 1. (Color online) Linear Bloch spectrum and Bloch waves at the symmetry points of the spectrum, (a) The linear spectrum 
for zero detuning. The Bloch waves in the lowest two bands are shown in (c) - (f), corresponding to the points indicated in 
(a), (b) The linear spectrum with a finite detuning. The corresponding Bloch waves are plotted in (g) - (j). Thick (red) lines 
are for (fii, while thin (blue) lines are for (f> 2 . The gray shadows represent — cos(2z) < 0. The unit of 2 : is 1/fciat. Note that we 
only show the lowest four bands. 


chemical potential. The amplitude spinor (j){z) satisfies III. LINEAR BLOCH SPECTRUM AND THE 
the time-independent GP equations PARITY SYMMETRY OF BLOCH WAVES 


, r 152 d S 1 , 

-^cos{2z)(l)++ \cj) 2 f)(l). (7) 

and one can see that Eq. (7) has a well-defined spin- 
dependent parity symmetry P if (5 = 0 


Since gap solitons only exist in the linear energy gaps, 
it is important to first identify the positions of these 
gaps. To do this, we use the Bloch ansatz (j){z) = 
Un,k{z) eyzp{ikz) in Eq. (7), while neglecting the non¬ 
linear terms. Here the Unk{z) are periodic functions, 
Un,k{z + tt) = Un,k{z)j u is the band index and k the 
quasi-momentum. The equations satisfied by the Un,k{z) 
are then 


P = P<J^, ( 8 ) 

where P is parity operator. This parity symmetry plays 
an important role in the specific form of the linear Bloch 
spectrum and therefore also for the properties of the gap 
solitons as we will show in the remainder of this paper. 

Once a gap soliton solution (j>{z) is found, its stability 
can be studied by simulating its long time behaviour us¬ 
ing Eq. (4). Eor this we add 10% Gaussian distributed 
noise to the wave function to trigger a possible instability 
and evolve the resulting state for 500ms. We also check 
the stability using the linear stability analysis by solving 
the standard Bogoliubov de Gennes equations. Both, the 
results from nonlinear evolution and the linear stability 
analysis, are consistent with each other. 




1 f d Y V 
( ^ + 1 Un,k - 2 COs{2z)Un,k (9) 
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dz 


ik 


s n 

+ -a. 


'^n,k 


and we calculate the linear Bloch spectrum by numerical 
diagonalization. 

The application of the parity symmetry operation 
given in Eq. (8) to above equations results in ^n,k = 
if (5 = 0. In fact, for these linear equations, the 
parity behavior is the same as the time-reversal-like sym¬ 
metry KcFx, with K being the complex conjugate opera¬ 
tor. While it might be natural to expect this symmetry in 
the spectrum, it is worth noting that it is broken by any 
non-zero detuning. This can also be clearly seen from the 
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FIG. 2. (Color online) Linear Bloch spectrum. Shaded re¬ 
gions represent the bands and white ones the gaps, (a) and 
(b): Spectrum as a function of the Rabi frequency fl for <5 = 0 
and <5 = O.SEiat, respectively. The blue vertical lines in (a) 
identify the values that will be used in Figs. 3 and 6, while the 
ones in (b) indicate parameters that relate to results shown 
in Figs. 8 and 9. (c) and (d): Spectrum as a function of the 
detuning S for fl = 8Eiat and fl — 2Eiat. The horizontal lines 
in (c) and (d) mark the values relating to results shown in 
Figs. 10 and 11, respectively. 

two typical spectra and associated Bloch wave solutions 
which are shown in Fig. 1: for zero detuning the spec¬ 
trum is symmetric [Fig. 1(a)], while for finite detuning 
an asymmetry is present [Fig. 1(b)]. 

From Figs. l(c)-(f) one can see that at the points 
around which the spectrum is symmetric (at fc = 0 and 
k = ±1), the Bloch waves for (5 = 0 possess the parity 
symmetry 


Vcjiiz) = X(l){z), (10) 

where the A = ±1 are the eigenvalues of the parity sym¬ 
metry operator V. This symmetry requires 

Re[(/)i(z)] = ARe[(() 2 (-^:)], Im[(()i(z)] = AIm[<() 2 (- 2 )]. 
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FIG. 3. (Color online) Gap solitons in the first and second 
gaps in the (N,^) plane for 5 = 0 and fl = 8Eiat. Shaded 
areas correspond to the linear bands. The soliton profiles at 
the labeled points in the first gap are shown in Fig. 4 and the 
ones in the second gap are shown in Fig. 5. 

If (5 < 0, the results for the first and second band are 
inverted. 

It is well-know that fundamental gap solitons can be 
understood as bifurcations from the linear Bloch waves 
at the symmetry points of the spectrum [3]. Consider¬ 
ing the intrinsic relationship between the fundamental 
gap solitons and the Bloch waves [4-7], we expect that 
the two features for the Bloch waves presented above are 
inherited by the fundamental solitons. 

The evolution of the linear Bloch band-gap structure 
as a function of the different parameters is given in Fig. 2 
by extracting the minima and maxima of ^n,k for each 
band in Fig.l(a)-(b). The shaded regions represent Bloch 
bands while the blank areas are the linear energy gaps. 
From the Fig. 2 the location of each gap can be easily 
determined. In following sections, we will study the exis¬ 
tence of gap solitons for the parameters indicated by the 
vertical and horizontal lines in Fig. 2. 


IV. GAP SOLITONS FOR THE CASE OF ZERO 
DETUNING 


It is very interesting to note that at these points the Bloch 
waves from neighboring bands have opposite parity. In 
Figs. I(c)-(f), the parity of the Bloch waves in the first 
band is A = —1, while the parity of the Bloch waves at 
the second band is A = 1. The sign of the parity of the 
Bloch waves depends on the sign of the Rabi frequency, 
fl. In the figures we show the solutions for O > 0, and 
for U < 0 the Bloch waves in the first band have A = 1 
and the second band A = — 1. If the detuning is finite 
(5 7 ^: 0), the Bloch waves become spin polarized and the 
polarization is determined by the sign of the detuning. 
Neighbouring bands have opposite spin polarization and 
in Fig. l(g)-(j), where we show results for 5 > 0, the 
second component (j )2 dominates in the first band while 
the first component (j)i dominates in the second band. 


In this section we study the gap solitons for the 
Rabi frequencies corresponding to the vertical lines in 
Fig. 2(a). The gap solitons are found numerically by 
solving Eq. (7) using the Newton relaxation method. One 
general feature of gap solitons in the absence of detuning 
{6 = 0) is that they are unpolarised {\(j)i{z)\'^ = \(j) 2 {z)\‘^) 
and solutions for the first and second gap for fl = 8Eiat 
are shown in Fig. 3. 

In the first gap only one family of fundamental modes, 
called the first family [5, 61], exists. Their spectrum is 
the lowest branch (dotted line) in the first gap and they 
are characterized by one main density peak confined in 
a unit cell as shown in Fig. 4(a). From Fig. 4(a) it can 
be seen that the first family has a parity symmetry with 
A = — 1. This family is expected to bifurcate from the 
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FIG. 4. (Color online) Soliton profiles at the points indicated 
in the first gap in Fig. 3. The grayed areas correspond to 
— cos(2z:) < 0. The thick (red) lines are for and the thin 
(blue) lines are for <j) 2 . 



FIG. 5. (Golor online) Time evolution of gap solitons in the 
first gap shown in Fig. 3. The initial wave functions are the 
gap soliton solutions with 10% Gaussian distributed noise 
added on. Plot (a) shows the out-phase mode of Fig. 4(c) 
and plot (b) shows the in-phase mode shown in Fig. 4(d). 
Only the density distribution of the hrst component, is 

shown. 

linear Bloch wave in the first band at the Brillouin zone 
edge [3] and it is very interesting to note that the parity of 
this family is same as the one of the corresponding linear 
Bloch wave (compare Fig. 1(d)). The nonlinear evolution 
demonstrates that the whole family is dynamically stable. 

Higher-order gap solitons in the hrst gap can be under- 
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FIG. 6. (Color online) Soliton profiles at the points indicated 
in the second gap in Fig. 3. (a) the hrst family of the funda¬ 
mental mode and (b) the second family. 

stood as composites from the fundamental mode. Three 
typical higher-order modes can be seen in Fig. 3 and the 
corresponding mode functions are shown in Fig. 4(b)-(d). 
The two dipole modes visible in Fig. 3 can be seen in 
Fig. 4(b) and (c) to correspond to a composite structure 
made from two fundamental modes in- and out-of-phase, 
respectively. The out-of-phase dipole modes are slightly 
more energetic than in-phase modes in the (N, fi) plane. 
The whole branch of the out-of-phase modes are unsta¬ 
ble, while the in-phase modes are stable except very near 
the linear Bloch bands. A typical unstable evolution of 
the out-of-phase mode triggered by the addition of an 
initial 10% Gaussian distributed noise is shown in Fig. 5. 
Additional higher-order modes can be constructed by in¬ 
cluding more the fundamental modes. The branch given 
by the dot-dash line in Fig. 3 in the first gap corresponds 
to solitons composed out of hve in-phase fundamental 
modes and its profiles are shown in Fig. 4(d). Note that 
the hve in-phase gap solitons grow in size and are stable 
except very close to the linear Bloch bands. This stable 
evolution is shown in Fig. 5. 

In the second gap we only consider the fundamental 
modes. There are two fundamental modes families, the 
hrst family (dotted line in Fig. 3 in the second gap) and 
the second family (solid line). Their prohles are shown in 
Fig. 6. The hrst family in the second gap can be consid¬ 
ered as a continuous extension of the hrst family in the 
hrst gap, as they have the same parity symmetry (i.e., 
A = — 1) and density distribution. The second family 
is particularly interesting, as they have double humped 
densities in the unit cell. This is very different from the 
case without spin-orbit coupling, where there is a node 
in density. Furthermore, the second family has a parity 
symmetry of A = 1, as can be seen from Fig. 6(b). The 
parity of the second family is the same as the one of the 
linear Bloch waves in the second band in Figs. l(e)-(f). 
The stability analysis demonstrates that the first family 
is stable except very close to the linear Bloch bands, while 
the second family is stable for most of the chemical po¬ 
tentials except a very narrow regime —3.9 < /i < —3.74 
where it is unstable. For comparison, we show an unsta¬ 
ble mode at fi = —3.75 and a stable mode /r = —3.7 in 
Fig. 7. 

For the case of H = 2Eiat, the first and second linear 
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FIG. 7. (Color online) Time evolution of the second family 
in the second gap as shown in Fig. 3. 10% Gaussian dis¬ 

tributed noise was added in the initial wave function to trigger 
potential instabilities. Plot (a) shows an unstable evolution 
a.t — —3.75 and plot (b) a stable evolution at /r = —3.7. 
Only the density distribution of the first component, \4>i\^, is 
shown. Note the different scale of t in (a) and (b). 



FIG. 8. (Color online) Fundamental families in the second 
gap for (5 = 0 and Q, = 2Eiat. Shaded areas represent the 
linear bands, (a) The dotted line corresponds to the first 
family and the solid line to the second family. The profiles 
of the labeled points are shown in (b)-(d). Since = 

\ 4 > 2 {z)'[^, only the densities of are shown and the thick (red) 
lines correspond to the first family, while the thin (blue) lines 
correspond to the second family. 

particle system in a harmonic trap. To do this, we start 
from corresponding Hamiltonian, 


bands overlap (see Fig. 2(a)). The first gap cannot open, 
but the second gap is large. The spectrum of two funda¬ 
mental families in this second gap is shown in Fig. 8(a) 
and, surprisingly, the two families are almost degener¬ 
ate in the (N, /r) plane and a crossover happens around 
/i = —2.75. The density distributions of these two fam¬ 
ilies (shown in Figs. 8(b)-(d)) are also very similar and 
only show small differences in the respective amplitudes. 
Both families are characterized by one main density peak 
within a unit cell. However, the parity symmetry of these 
two families takes the opposite sign just like in the sec¬ 
ond gap of Fig. 3(a), i.e., the first family has A = —I 
and the second family has A = 1. As in experiments the 
presence of matter-wave gap solitons is usually probed by 
observing the density distribution [12], the small density 
difference between these two families might not be suf¬ 
ficient to distinguish these two families with current ex¬ 
perimental abilities. For chemical potentials smaller than 
the crossover point, the second family is stable, while be¬ 
yond the crossover point it becomes unstable. The first 
family is stable except close to the third linear band. 

We have seen that for a large Rabi frequency (H = 
8 Eiat), the densities of the second family have a double 
humped structure, while for a smaller Rabi frequency 
(H = 2Eiat), the double humped structure disappears 
and the densities of the second family becomes single 
peaked. To understand this physically, we approximate 
the optical lattice potential of a unit cell by a harmonic 
oscillator potential and study in the following the energy 
spectrum and eigenstates of a spin-orbit-coupled single 


^^Har = + JPzCTz + —-f (II) 

Zm Z Z 

where is the harmonic trap frequency. We use the har¬ 
monic oscillator basis by replacing pz = iy^rnhjj^j2{a'^ — 
a) and z = y/h/2mu}z{o) + a) so that the Hamiltonian 
becomes. 


j. mhwz, is Al 

iiHar = nuJza'a + lyW —-— (o' - ajcTj, -I- —era,, (12) 


and we diagonalize it using hujz as the unit of energy. 
The ground and first excited states of this Hamiltonian 
are shown in Fig. 9, where one can clearly see that with 
increasing Rabi frequency the energy splitting between 
these two states increases. For a small Rabi frequency, 
the densities of both ground and first excited states have 
a single peaked structure [see Fig. 9(b)], but when in¬ 
creasing the Rabi frequency, the density peak of the first 
excited state becomes flat and gradually develops a dou¬ 
ble humped profile [see Fig. 9(c) and (d)]. Meanwhile the 
ground state keeps its single peak form. 

Approximating now a periodic optical lattice poten¬ 
tial by continuously connecting harmonic oscillator traps, 
the first and second families of the fundamental gap soli¬ 
tons in the lattice can be considered as extensions of the 
ground and first excited states in an individual harmonic 
trap respectively. Therefore some of the basic properties 
of the eigenstates of the harmonic trap will be passed on 
to the fundamental gap solitons and based on the profiles 
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FIG. 9. (Color online) Eigenspectrum and eigenstates of a 
spin-orbit-coupled single particle system in a harmonic trap 
with the Hamiltonian given in Eq. (12). (a) Lowest two 

eigenenergies. The dotted line is the ground state energy 
while the solid line is the energy of the first excited state. 
(b)-(d) Density distributions of corresponding eigenstates as 
indicated in (a). Since = 102 ( 2 ) 1 ^ we only show one 

component and the dotted lines are the densities of the ground 
state and solid lines are the ones o f the fir st excited state. 
Here 7 = 1.96 and the unit of z is . 

of the first excited states one can physically understand 
the appearance and disappearance of the double humped 
structure in the second family. Furthermore, we have 
also checked that the parity of the ground and first ex¬ 
cited states is same as that of the first and second families 
respectively. If we neglect the kinetic energy and without 
considering the optical lattice, our system in Eq. (7) be¬ 
comes a nonlinear Dirac equation [62] and we note that 
the double humped density is one of the features of soli- 
tons in the nonlinear Dirac equation [62]. 

Finally, we stress that the parity of the linear Bloch 
waves, the fundamental gap solitons and the eigenstates 
of the harmonic trap relates to the sign of the Rabi fre¬ 
quency. If the sign of the Rabi frequency changes, the 
parity switches sign accordingly. 


V. GAP SOLITONS FOR THE CASE OF 
NONZERO DETUNING 

The experimental preparation of a spin-orbit coupled 
BECs with zero detuning into the ground state is often 
accompanied by strong heating [13, 63], and finite de¬ 
tuning is known to help to suppress this effect. In this 
section, we study the case of nonzero detuning in which 
the spin composition of the gap solitons is unbalanced 
and parity symmetry is absent. To quantify the spin 
composition we define the spin polarization as 

a ^ [\Mz)\^ - \Mz)?] ,, 

Jdz[\Mz)\^ + \MzW]^ ^ ’ 

which is non-zero for d yf 0. 
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FIG. 10. (Color online) Fundamental gap solitons in the first 
and second gap for J = O.SEiat and D = 8 Eiat. (a) Soliton 
families in the (N, yj) plane. The dotted lines correspond to 
the first family and solid line to the second family. The shaded 
regions are the linear Bloch bands and in the embedded figures 
the spin polarization of corresponding families are shown. The 
soliton profiles at the labeled points are shown in (b) and (c), 
where the thick (red) lines are for (f>i and the thin (blue) lines 
are for ((> 2 . 

Gap solitons along the vertical line {6 = 0.5Eiat and 
D = SEiat) in Fig. 2(b) are shown in Fig. 10. Here we only 
show the fundamental modes, i.e., the first family in the 
first gap and the hrst and second families in the second 
gap. In the inset figures, the polarization of correspond¬ 
ing families is illustrated. The polarization of the first 
families in the first and second gap is positive while the 
polarization of the second family is negative. This can be 
confirmed by looking at the soliton profiles in Figs. 10(b) 
and (c). The second component ((^ 2 ) dominates in the 
first family (Fig. 10(b)) while the first component {(pi) is 
larger for the second family (Fig. 10(c)). The spin polar¬ 
ization of these two families coincides with that of linear 
Bloch waves in the first and second bands as shown in 
Fig. l(g)-(j). Soliton profiles of the first family in the 
first gap are same as the ones in Fig. 10(b). The second 
family still features a double humped density structure in 
each component, while the first family is single-peaked. 
The whole branch of the first family in the first gap is 
stable, while the first family in the second gap is stable 
in the regime ^ < —4.1, with the exception of points very 
close to the second linear band. Except for inside a nar¬ 
row regime given by —4.1 < /i < —3.6, the second family 
is stable. 

The gap solitons corresponding to the values indicated 
by the other vertical line (S = 0.5Eiat and O = 2Eiat) in 
Fig. 2(b) are shown in Fig. 11. A narrow first gap has 
opened due to the small detuning, which was previously 
closed in the linear spectrum for zero detuning. Since the 


































5=0.5E,, n=2E,, 
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n (units of EJ z (units of l/k, J 

FIG. 11. (Color online) Fundamental gap solitons in the sec¬ 
ond gap for S = O.SEiat and fl = 2Eiat. (a) Soliton families 
and (b) spin polarization. In both plots the dotted line corre¬ 
sponds to the first family and solid line to the second family. 
The shaded areas represent the linear Bloch bands. The soli¬ 
ton profiles of labeled points in (a) are shown in (c) and (d). 
The thick (red) lines are for </)i and the thin (blue) lines are 
for (/i 2 - The gray shadows indicate — cos(22;) < 0. 
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FIG. 12. (Color online) Behaviour of the Erst and second fun¬ 
damental families in the second gap with respect to detuning 
for fixed chemical potential = —4Ei, i.e., along the hori¬ 
zontal line in Fig. 2(c) with fl = 8E;. (a) Soliton families in 
the (N, S) plane and (b) spin polarization. In (a) and (b) the 
dotted lines correspond to the first family and the solid lines 
to the second family. The density distributions of the labeled 
points in (a) are shown in (c)-(h). The thick (red) lines are 
for (fii and the thin (blue) lines are for (j) 2 - The shaded areas 
represent — cos( 22 ) < 0. 

first gap is narrow, it is very difficult to find gap solitons 
there and we concentrate on studying the fundamental 
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n=2E,, ^=-3E,, 

lat ^ lat 



FIG. 13. (Color online) Behaviour of the first and second 
fundamental families in the second gap with respect to the 
detuning for fixed chemical potential jj, = —SE; along the 
horizontal line in Fig. 2(d). Here fl = 2Ei. 

modes in the second gap in which two families can be 
found. These have the same properties as the ones in the 
second gap of Fig. 10(a), however, compared to the case 
of zero detuning shown in Fig. 8, these two families do not 
overlap in the (N, /r) plane (see Fig .11(a)). Furthermore, 
the density distributions of the second family are different 
compared to the case of zero detuning: the component (j)i 
has a single main peak while the component (f )2 becomes 
double peaked. The stability calculations show that the 
whole second family is stable, while the first family is 
unstable for chemical potentials /r > —2.75. 

To further characterize the role of the detuning on the 
fundamental modes, we fix the chemical potential and 
study the existence of fundamental modes in the second 
gap with respect to a change in detuning for values indi¬ 
cated by the horizontal lines in Figs. 2(c) and (d). The 
results are shown in Figs. 12 and 13. For fl = 8Eiat 
[Fig. 12], the first family (dotted line in Fig. 12(a)) only 
exist for small detuning. With increasing values, the 
number of atoms in the first family increases, while it 
initially decreases and only later increases for the second 
family. The polarization of the second family also shows 
anomalous behavior for the large detuning regime. For 
large detuning the size of the second gap becomes very 
small [see Fig. 2], and the soliton profiles of the second 
family become less localized as shown in Figs. 12(g) and 
(h). For small detuning both families are stable, while 
they become unstable beyond the first extremum (around 
6 = lEiat) in their polarization curve. 

For fl = 2Eiat, as shown in Fig. 13, small values of the 
detuning can be seen to increase the spin polarization 
for the first and second families, while for large detuning 
both values saturate. Both family are stable except for 
detuning S > l.lEiat where the curvature of the polar¬ 
ization graph is changed. 


VI. CONCLUSION 

We have studied the properties of gap solitons in an 
experimentally available spin-orbit-coupled setup. When 
the detuning between the Raman beams and ground en- 
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ergy levels of the atoms is zero, parity symmetry plays 
important role. We have revealed that the two families 
of fundamental gap solitons in the second linear energy 
gap have opposite signs of parity, and the sign relates to 
the sign of Rabi frequency. The first family is character¬ 
ized by a single peaked density while the second family 
has a doubly humped density structure. However, the 
double humped structure disappears for small values of 
the Rabi frequency. An intuitive picture, that approxi¬ 
mates the lattice potential over a unit cell by a harmonic 
trap, was used to physically understand the evolution of 
the double humps. When the detuning is non-zero, fun¬ 
damental gap solitons become spin-polarized. The two 
families then have opposite polarization and the sign of 
polarization depends on the sign of the detuning. The 
parity symmetry and spin-polarization of these two fam¬ 
ilies of fundamental gap solitons correlate with that of 


linear Bloch waves at the symmetric points of the spec¬ 
trum. 

Finally, we would like to mention that beside the par¬ 
ity symmetry Eq. (7) also has a time-reversal-like sym¬ 
metry Kax if i5 = 0. If (j)o is a solution of this equation, 
then the wave-functions exp(i'!/i)(^o with arbitrary global 
phase 'i/' are also solutions and the uncertainty of the 
global phase makes a classification of solitons according 
to the time-reversal symmetry not straightforward [57]. 
We have therefore not discussed this case here. In ad¬ 
dition, Ref. [57] recently explored a further type of spin- 
dependent translational symmetry, which originates from 
spin-dependent optical lattices and which is also absent 
in our case. 

We acknowledge useful conversation with Avadh Sax- 
ena. This work was supported by Okinawa Institute of 
Science and Technology Graduate University. 
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